(Quasi)-exactly solvable quasinormal modes 
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We consider quasinormal modes with complex energies from the point of view of the theory of 
quasi-exactly solvable (QES) models. We demonstrate that it is possible to find new potentials 
which admit exactly solvable or QES quasinormal modes by suitable complexification of parameters 
defining the QES potentials. Particularly, we obtain one QES and four exactly solvable potentials 
out of the five one- dimensional QES systems based on the sl{2) algebra. 
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Introduction. Quasinormal modes (QNM) arise as per- 
turbations of stellar or black hole spacetimes They 
are solutions of the perturbation equations that are out- 
going to spatial infinity and the event horizon. Such 
"outgoing wave boundary condition" was first adopted 
by Gamow in his explanation of the a-decay of atoms as 
a quantum tunneling process . Generally, these condi- 
tions lead to a set of discrete complex eigenfrequencies, 
with the real part representing the actual frequency of os- 
cillation and the imaginary part representing the damp- 
ing. QNM carry information of black holes and neutron 
stars, and thus are of importance to gravitational-wave 
astronomy. In fact, these oscillations, produced mainly 
during the formation phase of the compact stellar ob- 
jects, can be strong enough to be detected by several 
large gravitational wave detectors under construction. 
Recently, QNM of particles with different spins in black 
hole spacetimes have also received much attention Q. 

Owing to the intrinsic complexity in solving the per- 
turbation equations in general relativity with the appro- 
priate boundary conditions, one has to resort to vari- 
ous approximation methods, eg., the WKB method, the 
phase- integral method etc., in obtaining QNM solutions. 
It is therefore helpful that one can get some insights from 
exact solutions in simple models, such as the inverted 
harmonic oscillator [1, 0] and the Poschl- Teller potential 
. Unfortunately, the number of exactly solvable models 
is rather limited. 

Recently, in non-relativistic quantum mechanics a new 
class of potentials which are intermediate to exactly solv- 
able ones and non-solvable ones has been found. These 
are called quasi-exactly solvable (QES) problems for 
which it is possible to determine analytically a part of the 
spectrum but not the whole spectrum 0, SH, [13) Ullj [ill ■ 
The discovery of this class of spectral problems has 
greatly enlarged the number of physical systems which 
we can study analytically. In the last few year, QES 
theory has also been extended to the Pauli and Dirac 
equations [13l] . 

In this paper we would like to study solutions of QNM 
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based on the Lie-algebraic approach of QES theory. We 
demonstrate that, by suitable complexification of some 
parameters of the generators of the sl{2) algebra while 
keeping the Hamiltonian Hermitian, we can indeed ob- 
tain potentials admitting exact or quasi-exact QNMs. 
Such consideration has not been attempted before in 
studies of QES theory. Our work represents a direct op- 
posite of the work in p^ . where QES real energies were 
obtained from a non-Hermitian PT-symmetric Hamilto- 
nian. 

QES Theory. Let us briefly review the essence of the 
Lie-algebraic approach 0, H, Q to QES models ■ Con- 
sider a Schrodinger equation Hip — E-tp with Hamiltonian 
H — —d^ + {dx = d/dx) and wave function ^p{x). 

Here x belongs either to the interval (— oo,cx3) or [0,oo). 
Now suppose we make an "imaginary gauge transforma- 
tion" on the function tp: ip{x) = x(x)e~^'^^-*, where g[x) is 
called the gauge function. For physical systems which we 
are interested in, the phase factor ex.p{—g{x)) is respon- 
sible for the asymptotic behaviors of the wave function 
so as to ensure normalizability. The function sat- 
isfies a Schrodinger equation with a gauge transformed 
Hamiltonian Hg = e^He^^. Suppose Hg can be writ- 
ten as a quadratic combination of the generators J" of 
some Lie algebra with a finite dimensional representa- 
tion. Within this finite dimensional Hilbert space the 
Hamiltonian Hg can be diagonalized, and therefore a fi- 
nite number of eigenstates are solvable. Then the system 
described by H is QES. For one-dimensional QES sys- 
tems the most general Lie algebra is sl{2), and Hg can 
be expressed as 

Hg^Y^ CabJ^J^ + XI + ^^^1 constant , (1) 

where Cab, Ca are taken to be real constants in [1, [^. 
The generators J" of the sl(2) Lie algebra take the dif- 
ferential forms: J+ = z^dz — nz , J° = zdz — n/2 , J~ — 
dz {n = 0,1,2,...). The variables x and z are re- 
lated by some function to be described later, n is the 
degree of the eigenfunctions x, which are polynomials 
in a (n -|- l)-dimensional Hilbert space with the basis 

(1,Z,Z2,...,Z")- 

Substituting the differential forms of J" into Eq. ([1]), 
one sees that every QES operator Hg can be written in 
the canonical form : Hg = —P4{z)dl + P3{z)dz + ^2(2), 
where Pk{z) are fc-th degree polynomial in z with real 
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coefficients related to the constants Cab and Ca- The 
relation between Hg and the standard Schrodinger oper- 
ator H fixes the required form of the gauge function g and 
the transformation between the variable x and z. Partic- 
ularly, X — j^' dy / ^Jl\ijj). Analysis on the inequivalent 
forms of real quartic polynomials P4 thus give a classifi- 
cation of all sZ(2)-based QES Hamiltonians [1, Q. If one 
imposes the requirement of non-periodic potentials, then 
there are only five inequivalent classes, which are called 
case 1 to 5 in 0. 

Our main observation is this. If some of the coefficients 
in Pk[z) are allowed to be complex while keeping V{x) 
real, then all the five cases classified in can indeed 
support QES/exact quasinormal modes. We shall discuss 
these cases below. 

QES QNM. We consider Case 3 in j9f, which cor- 
responds to Class I potential in Turbiner's scheme Q. 
There are two subclasses in this case, namely. Case (3a) 
and (3b). We shall present the analysis of QNM poten- 
tial for case (3a). The other case turns out to give the 
same potential with a suitable choice of the parameters. 
The potential in case (3a) has the form (in this paper we 
adopt the notation of @): 

V{x) = Ae^^^^ + Be^'= + Ce"^^ + De^^v^^ , (2) 

where x € (—00, 00) and v is a positive scale factor. Note 
that V{x) is defined up to a real constant, which we omit 
for simplicity, as it merely shifts the real part of the en- 
ergy. This remark also applies to the other cases. V{x) 
in Eq. ^ reduces to the exactly solvable Morse poten- 
tials when A = B = 0, or C = D = Q. This potential is 
QES when the coefficients are related by 



A 
C 



52 c+{n+l)iy I 
— , i> = 0, 

0,1,2... 



(3) 



Here 6, c, c? are arbitrary real constants. For each integer 
n > 0, there are n -I- 1 exactly solvable eigenfunctions in 
the [n + l)-dimcnsional QES subspace: 



i})n{x) = exp 



Xn{e 



2v 2^ 2v 

, ^ (4) 
Here Xn{z) is a polynomial of degree n in z = exp(Y^a;). 
To guarantee normalizability of the eigenfunctions, the 
real constants 6, c, d, u and n must satisfy certain rela- 
tions [i]. 

We want to see if we can get QNM solutions if we allow 
some parameters to be complex, while still keeping the 
potential V{x) real. This latter requirement restricts the 
possible values of the parameters, and hence the forms 
of QES potential admitting quasinormal modes. For the 
case at hand, we find that one possible choice of values 
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FIG. 1: Schematic sketch of the potential V{x) in Eq. (|6]) 
and the real parts of the corresponding QNMs with 6=1 and 
n = 1. 



of 6, c and d is: 

b = ib, c = — (n+ l)z/, d = d, b,d : real constants . (5) 
The potential Eq. ^ becomes 



Vn{x) 



[n + 1) de-^" + ^e-2V^^" , (6) 



and the wave function Eq. ([4]) becomes 



exp 



ib 



2v 

X Xn(e^") 



2v 



(7) 



V(x) approaches =Foo as x ^ ±00 respectively: it is un- 
bounded from below on the right. For small positive d 
and sufficiently large n, Vix) can have a local minimum 
and a local maximum. In this case, the well gets shal- 
lower as d increases at fixed value of n, or as n decreases 
at fixed d. Fig.[T]presents a schematic sketch of V(x) with 
b — \ and n = 1. We emphasize here that for different 
value of n, each V{x) represents a different QES poten- 
tial admitting n + 1 QES solutions. Since V{x) — > cx) as 
X —f —oo^ the wave function must vanish in this limit. 
This means d > from Eq. ([7]). For the outgoing bound- 
ary condition, we must take 6 > 0. Before we go on, 
we note here that there is another possible choice of the 
parameters, namely, 

b = —6, c ~ {n + d ^ id, b,d : real constants . (8) 

However, this choice leads to a potential related to Eq. ^ 
by the reflection x 1— > —x. Hence, we will only discuss 
the potential in Eq. ([6]) here. 

To see that the wave functions ipn{x) do represent 
quasinormal modes, we determine the corresponding en- 
ergy En- This is easily done by solving the eigenvalue 
problem of the polynomial part Xn{z) of the wave func- 
tion. From the Schrodinger equation we find that for 
n = 0, the energy is Eq = —v/A — ibd/2i'. This clearly 
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TABLE I: Values of QNM energy for the QES potentials in 
case 3 with parameters b = d = v = 1. Note that for each n, 
there are n+1 values of En- Energy levels of metastable states 
are in parentheses. WKB estimates of energy of metastable 
states are also listed. 



n 













-0.25 -0.5i 


— 


1 




(-2.349 - 0. 04491) 


-2.313 - 0.05881 






-0.151 - 0.955i 





2 




(-6.271 - 0.0001401) 


-6.267 - 0.0001471 






(-2.447 - 0.135i) 


-2.439 - 0.1561 






-0.0317 - 1.365i 




Q 
O 


\ 




1 9 9^^? 4 SI V 1 n~**i 






(-6.293 - 0.0008171) 


-6.284 - 0.0009021 






-2.542 - 0.2591 








0.0939 - 1.7401 




4 


(- 


-20.255 - 4.310 x lO'^^i) 


-20.254 - 4.328 x lO'^^i 
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(-6.323 - 0.0027421) 


-6.306 - 0.003111 
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shows that the only QES solution when n = is a QNM 
with an energy having a negative imaginary part (recall 
that b,d,i' > 0). For n = 1, we have two QES solutions. 
Their energies are Ei — —bv/A — ibd/2v ± \/v'^ — ibd. 
Again we have two QNM modes. One can proceed ac- 
cordingly to obtain n + 1 QNM modes with higher values 
of n. However, for large n, computation becomes tedious, 
and one has to resort to numerical means. For definite- 
ness, we list in Table H] some values of En for the case 
where h — d — v — 1. In this case, the potential has 
a local minimum and a local maximum, with the bar- 
rier height rises as n increases. For large n, some states 
(in parentheses) have their real parts of energy lie be- 
low the local maximum of the barrier . Such states are 
conventionally called the metastable states. We see from 
the table that the number of metastable states contained 
within the well increases as n becomes larger. This is 
reasonable as the well becomes deeper as n increases. 

We have also used the WKB method to determine pos- 
sible metastable states that can be trapped by the well. 
The method is as follows. Let the three turning points, 
from left to right, be designated as Xi,X2 and 2:3, as 
shown in Fig[T] The boundary conditions of QNM then 
leads to the following quantization condition for the en- 
ergy E: 



1 + 4e2T 



1 



where (3 
Qdx^) 



4e27 ' 
Qdx ^E-V{x), 



(9) 



and 



E. 



Re{E). This gives the quantization condition of Re{E): 
(3 = [l + l/2)7r , I = 0,1,2.... Then by keeping the 
imaginary part of in /5, but dropping it in 7, we can 
obtain an estimate of Im{E): 



Im{E) = - 



-2 j"^ dxy/V(x)-Re{E) 



2/;; dx/^R4E)-V{x) 



(10) 



In Table U we have also listed the WKB results for the 
metastable states. It is interesting to note that all the 
metastable states obtained by WKB methods are in fact 
the QES states in the cases we considered. The exact val- 
ues of the QES energies and those of WKB calculations 
are seen to be consistent. 

Oscillator-like potentials. We now turn to the other 
four cases. These four cases admit exact QNM solutions. 
Here we shall discuss Case 4 and 5 of [9|], which are as- 
sociated with oscillator potentials. 

For Case 4, the potential is given by 

V(x) = Ax*^ + Bx^ + Cx^ + ^ , xe\0,oo). (11) 

X'' 

It is QES if the coefficients are related by 
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with real constants b , c and d. As before, we now relax 
the reality constraint on the parameters but keeping V{x) 
real, and determine if QNMs can be supported in this 
case. It turns out the answer is positive, if we let b — 0, 
c = Aia, and d ~ d, with real a and d. This leads to the 
potential 



V{x) 



2 2 

-ax 



7(7-1) 



(13) 



To the first approximation, we set e '^'^ ~ and E 



Here 7 = d - § is arbitrary. Class VII and VIII in [8| 
also lead to this potential when some parameters are al- 
lowed to be complex. Eq. (fT^ is independent of n, and 
hence one can solve for solvable states with any degree n 
in its polynomial part. This system is therefore exactly 
solvable. For 7 > 1 and 7 < 0, V{x) is a monotonic de- 
creasing function on the positive half- line. If < 7 < 1, 
V{x) — > —00 at a; = 0, 00, and has a global maximum in 
between. 

The wave functions take the form = 
^ig^iax /2 ^^^^2^^ where Xn{x'^) is an n-th degree 
polynomial in z = x'^. In order that the wave function 
satisfies the outgoing wave condition and vanishes at the 
origin, we must have a > 0, 7 > 0. Plucking "ipnix) into 
the Schrodinger equation, we determine the energies to 
be En = — i(4n -|- 27 -|- l)a. Since a and 7 are positive. 
En is always negative, indicating decaying QNMs. 
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Case 3 and 4 illustrate the main steps in complexifying 
the relevant parameters to obtain QES/exact QNMs. It 
is satisfying to find that the remaining three cases can 
also be extended to give potentials which also admit exact 
QNMs. These cases will be treated briefly below. 

For Case 5 we find that the only viable choice of po- 
tential is V{x) = —{cx + d)^ + d^/4, where c and d are 
real constants. This is a shifted inverted oscillator. Just 
as in Case 4, here V{x) is also independent of n, and 
hence also exactly solvable. We mention that Class VI in 
Q also leads to this potential with appropriate choice of 
parameters. 

For simplicity, we briefly discuss the case with d = 0: 
V{x) — —0^x^/4: 4]. QNM solutions in such inverted 
oscillator has been discussed in Q using the modified 
annihilation and creation operators. Here we consider 
the problem from the point of QES potential. The wave 
function is given by ?A„(x) — exp(jcx^/4)x„(x). From the 
Schrodinger equation, we get the energies En — —icin + 
1/2). For c > 0, the wave function describes decaying 
outgoing QNM away from the maximum x = Q to x = 
±00. This is analogous to the QNM in black holes. When 
c < 0, we have growing incoming QNM moving towards 
the origin. This latter case was obtained in 

Hyperbolic potentials. Potentials in Case 1 and 2 in- 
volve hyperbolic functions. A proper complexification of 
the parameters in Case 1 is 

cd 

V{x) = — — tanh(-v/i/a;)sech(v^x) 

+ d^) sech2(y^a;) , (14) 

where x € (— (X),(X)). Its normalizable counterpart 
is, according to the classification in [l^, the exactly 
solvable Scarf II potential. The special case where 
d = 0, which is the inverted Poschl- Teller potential, 
has been employed in Q in their study of black hole's 
QNMs. We can easily obtain the wave functions and 
energies for the general case, Eq. (fT4|l . from QES 
theory. The wave function has the form: ipn{x) — 
(cosh 07a;)(*^+'')/2iv exp{id tan"! (sinh ^x) /2v) x 
Xn(sinh Y^a;), and the corresponding energy is 

^" = £-("+0 .-.c(n+i) . (15) 



Note that En is independent of d, which is a general fea- 
ture of Scarf- type potentials. As in Case 5, the imaginary 
part is proportional to n + 1/2, which is characteristic of 
black hole QNMs. We mention here that there is another 
complexification scheme for Case 1, giving a singular po- 
tential which exhibits very peculiar features, such as the 
existence of continuous bound states spectrum. As this 
is not related to QNM, it will be reported elsewhere [l7j. 

Case 2 potential with QNM is 



V{x) ~ — — coth(v^a;)cosech (v^a;) 

- —{v'^ + c^ + d^)coiiech^{^/T'x), (16) 

where x € (0,oo). Its normalizable coun- 
terpart is the generalized Poschl- Teller po- 
tential [lB|. The wave function is 4'n{x) — 
(sinh y/Dx) ^^'''rv)!2v (tanh(07a;/2))*''/2'^Xn (cosh ^x) . 
The energies are exactly given by Eq. (fTS]) . The two cases 
share the same QNM spectrum, despite the difference 
in the forms of the potentials and wave functions. This 
is not surprising, as it only reflects the same relation 
between the two counterpart potentials, the Scarf II and 
the generalized Poschl- Teller potential (see eg.. Table 1 
in [16]). 

To summarize, we have demonstrated that it is possible 
to extend the usual QES theory to accommodate QNM 
solutions, by complexifying certain parameters defining 
the QES potentials. We found that the five s^(2)-based 
QES systems listed in Q can be so extended. While one 
of these cases admits QES QNM, the other four cases 
give exact QNM solutions. It is hoped that our work 
would motivate the search of many more exact/quasi- 
exact systems of QNM in QES theories based on higher 
Lie algebras, and in higher dimensions. 
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